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ABSTRACT 

We present observations of ^^CO, C^^O, HCO+, Hi3C0+, DCO+ and N2H+ 
line emission towards the Barnard 68 starless core. The line profiles are inter- 
preted using a chemical network coupled with a radiative transfer code in order 
to reconstruct the radial velocity profile of the core. Our observations and mod- 
eling indicate the presence of complex radial motions, with the inward motions 
in the outer layers of the core but outward motions in the inner part, suggesting 
radial oscillations. The presence of such oscillation would imply that B68 is rel- 
atively old, typically one order of magnitude older than the age inferred from its 
chemical evolution and statistical core lifetimes. Our study demonstrates that 
chemistry can be used as a tool to constrain the radial velocity profiles of starless 
cores. 

Subject headings: astrochemistry - stars: formation — ISM: abundances — ISM: 
molecules — ISM: individual (Barnard 68) 



1. Introduction 



Starless cores likely represent the earliest stage of the formation of a star. In the standard 
view of star formation, cores form out an initially magnetically sub-critical cloud, and evolve 
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in a q uasi-static fashion through ambipolar diffusion ( Shu et aLlll987l : lMouschovias fc Morton 
199X1 ). In the opposite picture, cores are dynamic objects that form b y shocks in a turbu- 



lent flow, fragment and collaps e (or disappear) in a few crossing times (jPadoan et al.l 12001 



Ballesteros-Paredes et al.ll2003l ). Vel ocity measure r nents can provide important constraints 
on these two scenarios. For example, iTafalla et al.l (119981 ) examined the kinematics of L1544 
and inferred infall velocities up to 0.1 km s~^ a priori in compatible with sub-critical am- 
bipolar diffusion models. However, ICiolek fc Basul (120001 ) argued that these velocities can 
be understood if the core is super-critical. 

Line of sight velocities in cores can be established from observations of self-absorbed 
line profiles. In an infalling core, if the excitation temperature increases towards the center 
(as one can expect in a centrally condensed core with a roughly constant temperature), 
self-absor bed liiies are expected to be asymmetric, with a blue peak brighter than the red 
peak (see lEvanslll999l ). This techniqu e has been used in the past to detect collapse motions 
in starless cores and protostars (see iMyers et al.l l2000l . for a review). One of the major 
difficulties of this technique is the large degree of chemical complexity within these objects. 
For exam ple, numerous molecules are observed to deplete in the densest regions of starless 
cores (e.g. iBergin et al.ll2002l : iTafalla et al.ll2002l ). thus hampering our ability to measure the 
velocity in the dense central regions of these objects. 

Nevertheless, with a detailed knowledge of the chemistry of these objects, one can 
choose appropriate molecular tr ansitions to trace different part of th e cores, and unveil 
their radial kinematic structure (iBerginI l2003l : Ivan der Tak et al.l 120051 ) . In this paper, we 
use line observations of ^^CO, C^^O, HCO+, R^^CO^, DCO+ and N2H+ to reconstruct 
the velocity profil e of the Barnard 68 core (hereinaft er B68). With a well defined physical 



(Alves et al. 



2OOII: iBergin et al.l I2OO6I ) and chemical (IBergin et al.ll2002l : iMaret et al. 



2006 



Maret fc BerginI 120071 ) structure, this core is an ideal target for such a study. For this we 
compare the observed line profiles with the predictions of a chemical network coupled with 
a radiative transfer code. The paper is organized as follow. Observations are presented in 
[121 The analysis is detailed in ^ and the results are discussed in ^ 



2. Observations 

Maps of H13C0+ (1-0) {u = 86.754288 GHz), HCO+ (1-0) {u = 89.188523 GHz), 
N2H+ (1-0) {u = 93.173772 GHz) C^^O (1-0) {u = 109.782173 GHz), DCO+ (2-1) {u = 
144.077289 GHz), DCO+ (3-2) = 216.112604 GHz) and HCO+ (3-2) {u = 267.557526 
GHz) transitions were observed towards B68 {a = 17''22"^38.2'^ and 6 = -23°49'34.0"; J2000) 
in April 2002 and September 2002 with the IRAM-30m telescope. The ^^CO (2-1) {1^ = 



- 3 - 



220.398684 GHz) transition was observed in September 2003 with the Caltech Submi l limete r 
Observatory (CSO) . These observations have aheady been presented in iBergin et al.l (120021 ) , 
Bergin et al.l (120061 ). iMaret et al.l (120061 ) . and lMaret &: BerginI (120071 ). and we refer the reader 
to these papers for technical details. 

Transitions of HCO+ (4-3) {u = 356.734134 GHz) and N2H+ (3-2) {u = 279.511832 
GHz) were observed in August 2005 using the Atacama Pathfinder eXperiment telescope 
(APEX). The APEX-2A receiver was used together with the FFTS spectrometer, which 
provides a spectral resolution of 61 kHz. The telescope half power beam size is 22" at 280 
GHz and 17" at 357 GHz. The data were calibrated using the chopper wheel method, and the 
system temperature was ~100 K at 280 GHz and 130-180 K at 357 GHz. The observations 
were converted on the scale assuming a main beam efficiency of 70% and a forward 
efficiency of 95% (C. de Breuck, priv. comm.). These observations were obtained in beam 
switching mode. 

Fig. [1] present the observations together with our model predictions (see § |3] for a 
description of the model). The lines show a variety of different profiles, from Gaussian to 
self-absorbed and asymmetric double peaked profiles. C^^O, ^^CO and H^'^CO+ are probably 
optically thin (or only moderately optically thick), and have centrally peaked Gaussian 
profiles. On the other hand, HCO+ (1-0) and (3-2) lines are self-reversed and show the 
typical blue asymmetry (blue peak brighter than the red peak). These line profiles suggests 
that collapsing motions are present along the line of sight. Surprisingly, self-absorbed HCO'*" 
(4-3), DCO+ (3-2), and the central component of the N2H+ (1-0) line exhibit red asymmetry 
(red peak brighter than the blue peak), indicating expansion. Finally, the DCO^ (2-1) line 
is self-absorbed but symmetric, hinting at the presence of static gas along the line of sight. 
The absorption dip velocity of self-reversed lines is consistent with the wlsr of the source 
(3.31 kms~^), as determined from the C^^O (1-0) line. In the following, we model the line 
profiles in order to constrain the line-of-sight velocity of the core. 



3. Analysis and results 



We h ave modelled t he sp e ctral line profiles sh own on Fig. [T] following the approach de- 
scribed in iBergin et all kood \. iMaret et all kood \ and iMaret fc BerginI J2007h . The abun- 
dance profiles were computed using a chemical network that includes gas grain interactions, 
as we ll as the fractionation reactions for carbon, oxygen, and deuterium (IBergin fc Langer 



19971 ). These abundance prof iles were used to c ompute the line emission using a 1-D Monte 
Carlo radiative transfer code (lAshby et al.ll2000l ). These model predictions were convolved at 
the appropriate spatial resolution in order to be compared with the observations, assuming 



-4- 



a distance of 125 pc. 



The den s ity p rofile from lAlves et al.l (120011 ) and the gas and dust temperature from 
Bergin et al.l (120061 ) were used. We have adopt ed the same initial abundances and cosmic 
ionization rate as in the iMaret fc Bergiru (120071 ) "best-fit" model. As a first approach, we 
have assumed that no systematic radial motions are present in the core. We have used the 



turbulent velocity profile determined by iBergin et al. 



mom from C^^O (1-0) fine observa- 



tions. The turbulent velocity contribution is 0.3 km s~^ (FWHM) at the surface of the core, 
and decreases towards the center down to 0.15 km s~^. Rotation of the cloud has been 
neglected. Both the ^^CO (2-1) and the H^'^CO^ (1-0) lines have an unresolved hyperfine 
structure that need to be taken into account i n ord er to properly model the line profiles. 
Following the approach used by iTafalla et al.l (120061 ) , we have artificially broadened these 
two lines by 0.13 km s^^, w hich corresponds to the s eparation between the two hyperfine 
components of each species (ISchmid-Burgk et al.ll2004l ). 



On Fig. [21 we present the abundances predicted by our chemical model. From this figure 
we can roughly estimate which lines originate in different parts of the core. For example, the 
HCO"*" abundance peaks in the outer part of the core {A^ ~ 2) and decreases towards the 
center because of the freeze-out of its parent molecule CO. Consequently, HCO+ lines probe 
predominantly the outer layer. Both N2H"'' and DCO'^ abundances peak deeper inside, at 
an Ay of 5 and 13 respectively. Therefore lines from these two species preferentially trace 
the inner regions of the core. Of course, excitation effects are also important, and different 
transitions of the same molecule are expected to probe different depths in the core as well. 
For example, the HCO"'" (4-3) opacity is likely lower that the (1-0) and (3-2) lines, and thus 
the HCO"'' (4-3) line probably arises from deeper regions than the HCO"*" (3-2) and (1-0) 
lines. Altogether, our observations suggest the presence of complex radial motions in B68, 
with inwards motions in the outer part of the core - as indicated by the blue HCO^ (1-0) 
and (3-2) line profiles - and outward motions in the inner region - as shown by the HCO"'" 
(4-3), DCO+ (3-2) and N2H+ (1-0) lines. In addition, the DCO+ (2-1) line hints at the 
presence of an intermediate static region between these two parts. 

In order to test this hypothesis, we have computed the line profiles for the step-like 
velocity profile represented on Fig. [31 By convention, negative velocities correspond to 
inward motions, while positive velocities correspond to outward motions. The velocity is 
assumed to be negative (-0.045 km s~^) at radius greater than 8,000 AU, and positive (0.025 
km s~^) at smaller radii. Fig. [1] compares the observations with our model predictions for 
this velocity profile. The fit is not perfect. However, given the complexity of our modeling - 
chemical network coupled with a radiative transfer model - the overall agreement between 
the model predictions and observations is fairly good. The model reproduces quite well the 
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infall asymmetry of HCO"'" (1-0) and (3-2), as well as the blue asymmetry of the (4-3) line. 

(1-0), N2H+ (3-2), C^^O (1-0) and ^^CO (2-1) hne profiles are also well matched. 
Although the model predicts correct integrated intensities for DCO"'" (3-2) and N2H+ (1-0), it 
does not reproduce the blue asymmetries seen in these lines. The predicted line profiles have 
a Gaussian shape, whereas the observed lines are self-absorbed. In these cases, it appears 
that our model underestimates the opacity of these lines. 



4. Discussion and conclusions 



Our observations and modelling indicate that the outer part of the core is collapsing 
while the inner part is expanding. The velocities in both regions are relatively small (a few 
tens of meters per second) and are largely sub-sonic. The transition between the collapsing 
and expanding region occurs at a radius of about 8,000 AU. It is important to note that the 
fit to the observations is probably not unique. However, large changes in the transition and 
velocities between the regions are likely ruled out. As already mentioned, our model does 
not reproduce the asymmetries of the N2II+ (1-0) and DCO+ line profiles, due to an under- 
estimating of the opacity. There are several possible explanations for this discrepancy. First, 
our chemical model might not predict the correct abundance profile for these two species. 
For example, increasing the DCO"*" abundance at intermediate radii, with a subsequent re- 
duction at small radii to keep the column density constant would increase the opacity of 
line without changing the line flux significantly (which is correctly predicted). Indeed, our 
model predic ts a DCO''' (2-1) line intensity at 10 < < 20 that is lower than observed 



(see Fig. 3 in lMaret fc Berginll2007l ). Second, the physical structure of the core is also some- 
what uncert ain. For exaniple, th e ter nperature profile rn ay be slightly different than the one 
predicted by iBergin et al.l (120061 ) (see lCrapsi et al.l 120071 ). although small changes (±3 K) in 
the temperature profile were found to have little effect on the line intensities. Turbulence at 
the center of the c ore might also be di fferent that what assumed here. The turbulent profile 
was determined by lBergin et al.l (120061 ) using C^^O observations, which are not a good probe 
of the innermost regions of the core because of heavy depletion. Decreasing the turbulent 
velocity profile in the innermost region of the core would also increase the opacity of these 
two lines. Regardless of the discrepancy with our model predictions, we emphasize these two 
line profiles, together with the HCO^ (4-3) line profile, unambiguously indicate the presence 
of outwards motions at the center of the core. We found that no constant infalling velocity 
profile could reproduce the observations. 



The step-like velocity profile we obtain is suggestive of the presence of radial oscillations 
in the core. Small non-radial oscillations of the outer layers of the core around a stable equi- 
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librium have been proposed by iLada et al.l (120031 ) and iRedman et al.l (120061) to interp r et the 
asymmetry of the CS (2-1) and HCO"'" (3-2) hnes observed in Barnard 68. iLada et al.l (120031 ) 
also suggested that, in addition to the surface oscillations seen in the CS (2-1) line, some 
radial oscillation might also be present. Our observations and modeling are in agreement 
with this hypothesis. 



Keto fc Fieldl (120051 ) proposed a hydrodynamic model for the evolution of cores in quasi- 
equilibrium. In their modeling, a small perturbation of a stable core (according to the 
Bonnor-E bert criteri o n) is found to result in damped oscillations around the equilibrium 
position. iKeto et al.l (120061 ) showed that the surface pattern observed in the CS (2-1) line 
towards B68 is consistent with the (n = 1, 1 = m = 2) quadripolar oscillation mode for a 
specific orientation {X = 9 = 30°). For the same oscillation mode, the radial velocity is 
qualitatively similar to the one obtained in this study, but is out of phase: the inner part of 
the core has a negative radial velocity, while the outer part of the core has a positive velocity. 
At the other half of the oscillation cycle, the velocity would be reversed and in qualitative 
agreement with the profile we obtain (E. Keto, priv. comm.) 



As pointed out by iKeto et al.l (120061 ) and iRedman et al.l (120061 ) . the presence of oscilla- 
tions in B68 suggests that the core is relatively old, typically older than a few sound speed 
crossing times (~ 10^ yr). This is comparable to the ambipolar diffusion timescale at the 
center of the core and more tha n one order of magnitu de higher than the free-fall time scale 
(10^ and 7 x lO'* yr respectively; iMaret fc Berginll2007l ). This is also about an order of mag- 



nitude higher than the ag e of the core as determ ined from the degree of CO depletion (lO^yr; 



Berlin et al. 


2002. 


2006 





For t = 3 X 10^ 



Bergin et al.l (120061 ) model predicts a C^^O 
(1-0) line intensity that is three time smaller than the observations, while at t = 10^ yr, the 
greater CO depletion leads to an even greater mismatch between observations and model 
predictions. The N2H+, HCO+ H^^CO^, and DCO+ would also be underestimated by the 
model because of the freeze-out of their parent of their precursors, N2, CO and ^'^CO. Thus 
the "chemical age" of B68 inferred from the CO depletion appears to be hard to reconcile 
with the "dynamical age" implied by the presence of surface oscillations. It is also interesting 
to compare the sound crossi ng time in B68 to typical statistical starless core lifetimes (see 
Ward-Thonapson et al.ll2007l . for a review). Using sub-millimeter continuum emission maps. 



Kirk et al.l (l2005l ) obtained lifetimes of 1 — 3 x 10^ years which are in good agreement with 
the age of the B68 determined from CO depletion, but again about an order of magnitude 
higher than the sound crossing time. 



^Strictly speaking, the age determined by this method is a lower hmit of the real age of the clou d, because 
CO is assumed to be pre-existing, and the density is supposed to be constant (|Bergin et all 120061 ). 
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The present study emphasizes the importance of understanding the chemistry to con- 
strain the kinematics of starless using spectral line profiles. Because of the strong chemical 
gradients that exist in these cores - as well as excitation effects - different lines can be used 
to selectively trace different parts of the core. We have demonstrated that, using lines that 
are appropriately chosen, it is possible to reconstruct the core radial velocity profile. In B68, 
our observations and modelling suggest the presence of complex line-of-sight motions that 
are consistent radial oscillations. The presence of such oscillation indicates that some core 
are long-lived, with lifetimes about an order of magnitude higher than those derived from 
their chemical evolution or sub-mill imeter surveys. S o far, oscillations have been inferred 
in only one other core, FeSt 1-457 (lAguti et al.l 120071 ). Clearly, similar studies on a larger 
number of cores are needed to establish if such cores are common. 



S. M. wishes to thank Eric Keto for fruitful discussions. The authors are also grateful 
to the referee for useful and constructive comments. This work is supported by the National 
Science Foundation under grant 0335207. 
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Fig. 1. — Comparison between the line profiles observed towards the extinction peak of B68 
(black histograms) and the predictions of our "best fit" model (sohd grey hnes). For the 
clarity of the plot, lines are shifted vertically, and some of them are scaled by the factor 
mentioned on the right of each curve. 
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Fig. 2. — Abundances predicted by our chemical model for the observed species as a function 
of the visual extinction in the core. 



- 12 - 



0.04 



0.5 



0.02 



-0.02 



-0.04 







5,000 10,000 
Radius (AU) 



- 0.4 



0.3 



0.2 "'a 



0.1 



Fig. 3. — Best fit radial velocity profile. The solid line shows the systemic velocity profile 
(values are on the left axis). The dashed l ine shows the turbulent velocity profile (values are 
on the right axis, from lBergin et al.ll2006l ). 



